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Abstract 

An elimination problem in semidefinite programming is solved by means of 
tensor algebra. It concerns families of matrix cube problems whose constraints 
are the minimum and maximum eigenvalue function on an affine space of sym- 
metric matrices. An LMI representation is given for the convex set of all feasible 
instances, and its boundary is studied from the perspective of algebraic geometry. 
This generalizes the earlier work [12] with Parrilo on fc-ellipses and fe-ellipsoids. 
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1 Introduction 

The matrix cube problem in semidefinite programming is concerned with the 
following question. Given real symmetric x A^-matrices Aq, Ai, . . . , Am , does 
every point (ti, . . . , tm) in the cube YYlLi /"i] satisfy the matrix inequality 

m 

Ao + Y,tkAk h ? 

k=l 

The inequality means that the symmetric matrix Aq + YlT=i ^k^k is positive 
semidefinite, i.e. its N eigenvalues are all non-negative reals. This problem is NP- 
hard [3], and it has important applications in robust optimization and control 
[5 [H]) e.g., in Lyapunov stability analysis for uncertain dynamical systems, 
and for various combinatorial problems which can be reduced to maximizing a 
positive definite quadratic form over the unit cube. For a recent study see [5]. 

In this paper we examine parametrized families of matrix cube problems, 
where the lower and upper bounds that specify the cube are the eigenvalues of 
symmetric matrices that range over a linear space of matrices. We define the set 



C 



(x, d) G X 



m 

d ■ Aq -\- ^ tfc^fe ^ whenever 
1 — 1 



k=l 

.mm{Bk{x)) < tk < Amax(-Bfc(2;)) 

for k = 1,2, ... ,m 



(1.1) 
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where the Ai are constant symmetric matrices of size Nq x Nq, the symbols 
Amin(') and Amax(') denote the minimum and maximum eigenvalues of a matrix, 
and the x Nk matrices Bk{x) are linear matrix polynomials of the form 

Here B^^\ b[^\ . . . , B^it'^ are constant symmetric N}^ x Nj^ matrices for all k. 

Since the minimum eigenvalue function Amin(') is concave and the maximum 
eigenvalue function Amax(') is convex, we see that C is a convex subset in M""*"^. 
Our definition of C is by a polynomial system in free variables (x, d) and uni- 
versally quantified variables t. Quantifier elimination in real algebraic geometry 
[l] tells us that C is semialgebraic, which means that it can be described by a 
Boolean combination of polynomial equalities or inequalities in {x,d). However, 
to compute such a description by algebraic elimination algorithms is infeasible. 

Linear matrix inequalities (LMI) are a useful and efficient tool in system and 
control theory [2]. LMI representations are convenient for building convex opti- 
mization models, especially in semidefinite programming [16] . The set described 
by an LMI is always convex and semialgebraic. Since our set C is convex and 
semialgebraic, it is natural to ask whether C admits an LMI representation ? 

Our aim is to answer this question affirmatively. Theorem 13.21 states that 

C = {(x,d) € M" X M : /:(x,d) ^ 0}, (1.2) 
where C{x,d) is a linear matrix polynomial whose coefficients are larger symmet- 

(k) 

ric matrices that are constructed from the matrices Ai and Bj . The construc- 
tion involves operations from tensor algebra and is carried out in Section 3. 

First, however, some motivation is needed. In Section 2 we shall explain why 
the convex semialgebraic sets C are interesting. This is done by discussing several 
geometric applications, notably the study of m-ellipses and m-ellipsoids [T2] . 

Section 4 is devoted to algebraic geometry questions related to our set C: 
What is the Zariski closure Z{dC) of the boundary dC ? What is the polynomial 
defining Z{dC) ? What is the degree of Z{dC) ? Is the boundary dC irreducible ? 

In Section 5, we discuss an application to robust control, explore the notion 
of matrix ellipsoids, and conclude with some directions for further research. 



2 Ellipses and beyond 

In this section we illustrate the construction of the set C for some special cases. 
The first case to consider is A^i = • • • = = 1 when each Bi^{x) = b'j^x + (3^ is 
a linear scalar function. Then our elimination problem is solved as follows: 

C = |(x,(i)GM"+i: (i-Ao + f](6^x + /3fc)Afcho|. 

Thus C is a spectrahedron (the solution set of an LMI) and, conversely every 
spectrahedron arises in this way. The Zariski closure Z{dC) of its boundary dC is 
the hypersurface defined by the vanishing of the determinant of the above matrix. 
For generic data Ai,bk,Pk this hypersurface is irreducible and has degree A^o- 
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Throughout this paper, we shall use the term "generic" in the sense of algebraic 
geometry (cf. |13j). Randomly chosen data are generic with probability one. 

A natural extension of the previous example is the case when each given 
matrix Bk{x) is diagonal. We write this diagonal matrix as follows: 

B,{x) = diag(6i')^x + (3l^\ h^^'^^^x + (3['''^) 
Then our spectrahedron can be described by an intersection of LMIs: 

The Zariski closure Z{dC) of the boundary dC is a hypersurface which is typically 
reducible. It is defined by the product of all determinants of the above matrices 
which contribute an active constraint for C. Each of the determinants has degree 
A'^o and there can be as many as Ai A2 • • • of these boundary components. 

The point of departure for this project was our paper with Parrilo [12] whose 
result we briefly review. Given points (ui, fi), . . . , {um,Vm) in the plane and 
a parameter d > 0, the corresponding m-ellipse £d is the set of all points (xi, X2) 
whose sum of distances to the given points (uj, Vi) is at most d. In symbols, 

£d = |x G R2 . ^(xi - UkY + {X2 - VkY < dj . 

In |12j it is shown that 8^ is a spectrahedron, an explicit LMI representation of 
size 2"^ X 2"* is given, and the degree of d£ is shown to be 2™ when m is odd, and 
2™ — (^2) when m is even. For instance, if m = 3, (ni,fi) = (0,0), (1*2,^2) = 
(1,0) and (lis, V3) = (0, 1), then the 3-ellipse Ed has the LMI representation 



■f 3x1 - 1 


X2 - 1 


X2 





X2 











X2 - 1 


i + xi - 1 





X2 





X2 








X2 





d + xx + l 


X2 — 1 








X2 








X2 


X2 — 1 


d — xi + I 











X2 


X2 











d + xi -1 


X2 — 1 


X2 








X2 








X2 — 1 


d — xi — 1 





X2 








X2 





X2 





d — xi + 1 


X2 — 1 











X2 





X2 


X2 — 1 d 


- 3X1 + 



>- 0. 



See Figured] for a picture of this 3-ellipse and the Zariski closure of its boundary. 

We can model the m-ellipse using matrix cubes as in (|l.ip as follows. Let 
Nq = 1 and define each scalar Ai to be 1, and let A'^i = • • • = Nm = 2 and define 



Bk{x) 



Xl -Uk X2- Vk 
X2 - Vk Uk- Xl 



for /c = 1, • • • , m. 



The eigenvalues of Bk{x) are ± {xi — UkY + (2^2 — Vk)"^-, and we find that 

C = {{x,d) & : d>ti + --- + tm whenever \tk\ < ^ {xi-UkY + {x2-VkY}. 

This formula characterizes the parametric m-ellipse 

C = {{xi,X2,d) ■ {xi,X2)e£d], 
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Figure 1: The Zariski closure of the 3-elhpse is a curve of degree eight. 



and its LMI representation (|1.2p is precisely that given in |121 Theorem 2.3]. The 
construction extends to m-ellipsoids, where the points lie in a higher-dimensional 
space. We shall return to this topic and its algebraic subtleties in Example 14.41 
In Section 5 we introduce a matrix version of the m-ellipses and m-ellipsoids. 

Consider now the case when Nq = 1 and Aq = Ai = ■ ■ ■ = Am = — 1 but the 
-Bfc(x) are allowed to be arbitrary symmetric A^^^ x A'^fc-matrices whose entries are 
linear in x. Then the spectrahedron in (jl.ip has the form 

m 

C = gM" xM : ^Amax(Sfc(2;)) <d}. 

k=l 

A large class of important convex functions on can be represented in the 
form X I— > Amax(-Bfc(2;)). Our main result in the next section gives a recipe for 
constructing an explicit LMI representation for the graph of a sum of such convex 
functions. The existence of such a representation is not obvious, given that the 
Minkowski sums of two spectrahedra is generally not a spectrahedron [141 §3.1]. 

3 Derivation of the LMI representation 

In this section we show that the convex set C in (jl.ip is a spectrahedron, and we 
apply tensor operations as in [12] to derive an explicit LMI representation. We 
begin with a short review of tensor sum and tensor product of square matrices. 
Given two square matrices F = {Fij)i<ij<:r and G = {Gki)i<k,e<s, their standard 
tensor product ^ is the block matrix of format rs x rs which defined as 

F (g) G = {FijG)l<ij<r- 

Based on tensor product 0, we define the tensor sum © to be the rs x rs-matrix 

F®G = F0ls + It(8)G. 
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Here denotes the identity matrix of size r x r. For instance, the tensor sum 
of two 2 X 2-matrices is given by the formula 

"a + e / 6 



a b 




"e / 




g a + h Ob 


c d 




g h 




c d+e f 



c g d + h 

Tensor products and tensor sums of matrices are also known as Kronecker prod- 
ucts and Kronecker sums [9] . Note that the tensor operations (8> and © are not 
commutative but they are associative. Hence we can remove parentheses when 
we take the tensor product or tensor sum of k matrices in a fixed order. 

The eigenvalues of F © G are the products of pairs of eigenvalues of F and G. 
Similarly, the eigenvalues of FffiG are the sums of such pairs. This is well-known 
for tensor products, but perhaps slightly less so for tensor sums. We therefore 
explicitly state the following lemma on diagonalization of tensor sums. 

Lemma 3.1. [12^ Lemma 2.2] Let Mi, . . . , Mf^ be symmetric matrices, Ui, . . . ,Uk 
orthogonal matrices, and Ai, . . . , diagonal matrices such that Mi = Ui-Ki-Uf 
for i = 1,. . . ,k. Then the tensor sum transforms as follows: 

(C/i © ■ ■ ■ (g) Uk)^ • (Ml © ■ ■ ■ © Mfc) ■ (?7i © • • • © t/fc) = Ai © • • • © Afc. 

In particular, the eigenvalues o/Mi©M2©- • -©Mfc are the sums Ai + A2 + - • •+Afc 
where Ai is any eigenvalue of Mi, A2 is any eigenvalue 0/M2, etc. 

We now turn to the construction of an LMI representation for the convex 
semialgebraic set C. First, we introduce a linear operator A on the space of 
linear matrix polynomials in m + 1 unknowns a = (oq, ai, . . . , flm)- Namely, let 
Pa{x,d) be such a linear matrix polynomial of the form 

Pa{x,d) = aQPQ{x,d) + aiPi{x,d) ^ VamPm{x,d) (3.1) 

where the coefficients PQ[x,d), . . . , Pm{x, d) are matrix polynomials in (x, d). Us- 
ing the matrices Aq, Ai, . . . , An in (jl.ip . we define the linear operator A by 

A{Paix,d)) := Poix,d)®Ao + Pi{x,d)(g)Ai + --- + P„,{x,d)<^Am. 

Second, using the matrices Bk{x) in (II. 1|) we define the following tensor sum 

La{x,d) := {dao)e{aiBi{x))e ■■■ ®{amBmix)). 

This expression is linear matrix polynomial in a = (ao, ai, . . . , Om), and we set 

C{x,d) := A{La{x,d)). (3.2) 

Since La{x, d) is linear in both (x, d) and a, and since .A is a linear operator, the 
matrix £(x, d) depends linearly on (x, d) and on the matrices ^i) • • • > ^m- 

Theorem 3.2. The convex semialgebraic set C C M*^"*""^ in lll.l]) is a spectrahe- 
dron, and it can be represented by the LMI il.2\) . 
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This theorem imphes that optimizing a hnear function over the convex set C 
is an instance of semidefinite programming |16j . For instance, for a fixed value 
d* of d, to minimize a linear functional c^x over its section C H { (x, d) : (i = d*} 
is equivalent to solving the semidefinite programming problem 

minimize (Fx subject to £{x,d*) ^ 0. 

This problem can be solved by standard SDP solvers like SeDuMi |15j . 

To prove Theorem 13.21 we need one more fact concerning the linear operator 
A, namely, that A is invariant under congruence transformations. 

Lemma 3.3. Let Pa{x,d) be as in 13. For any matrix U{x,d) we have 

A{U{x,dfPaix,d)Uix,d)) = {U{x,d)0lNofA{Paix,d))iU{x,d)^lNo)- (3.3) 

Proof. First note that tensor product satisfies {A <S> B)'^ = and 

(Ml M2) • (M3 ® M4) • (M5 ® Me) = (Ml • M3 • M5) (g) (M2 • M4 • Me) 
Using these identifies we perform the following direct calculation: 

{u{x,d)®lNoY A[Pa{x,d)^ (uix,d)0lN^^ 
= (u{x, d) ® Ino) Y1 d-)®^k\ (U{x, d) ® Ino) 

\k=0 J 

= J2(uix,d)(S)lNo) (^Pkix,d)®Ak^(u{x,d)®lNo) 

k=0 

m 

= Y,{u{x,dfPk{x,d)U{x,d)^ ®Ak. 

k=0 

By definition of A, this expression is equal to the left hand side of (j3.3p . □ 
Proof of Theorem \3.2l Let Uk{x) be orthogonal matrices such that 

Ukixf ■ Bk{x) ■ Uk{x) = diag(Af)(x),...,ASJ^^(x)) =: Dfc. (3.4) 

Here X^^\x) are the algebraic functions representing eigenvalues of the matrices 
Bk{k). If we set Q = (l) ® Ui{x) ® ■ ■ ■ ®i Um{x) then, by Lemma [3TH we have 

■ La{x,d) -Q = {dan) ® {aiDi{x)) ® ■ ■ ■ ® {araDm{x)) =: L. 

Note that L is a diagonal matrix, with each diagonal entry having the form 

m 

d-ao + y^^xfj{x)ak. 

k=l 

It follows that A{L) is a block diagonal matrix, with each block of the form 

m 

d-Ao + ^xfJ{x)Ak. (3.5) 

k=l 
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By Lemma 13.31 and the definition ol C{x,d), we have 



A(q'^ ■ La{x,d) -Q) = A{L). (3.6) 



I No ) 



Hence C{x,d) ^ and only if the N1N2 ■ ■ ■ Nm blocks (|3.5p are simultaneously 
positive semidefinite for all index sequences ji,j2, ■ ■ ■ ,jm- Given that the ^^Jl^i^) 
are the eigenvalues of Bk{x) we conclude that C{x, f^) ^ if and only if 



d ■ Aq + tkAk ^ for all (ti, . . . , t^) such that is eigenvalue of -Bfc(a 



k=l 



Since d ■ Aq + XlfcLi ^kAk h describes a convex set in t-space, we can now 
replace the condition "t^, is an eigenvalue of Bk{x)" by the equivalent condition 



^mm{Bk{x)) < tk < Amax(-Bfc(2;)) for /c = 1, 2, 

We conclude that C{x,d) ^0 if and only if (x, d) S C. 



, m. 



□ 



The transformation from the given symmetric matrices Ai and Bj(x) to the 
bigger matrix C{x,d) is easy to implement in Matlab, Maple or Mathematica. 
Here is a small explicit numerical example which illustrates this transformation. 



Example 3.4. Let n = m = Nq = Ni 

Ao 



"2 


l' 


, Ai = 


"1 


l' 


, A2 = 


"0 


1" 


1 


2 


1 





1 


1 



N2 = 2 and consider the input matrices 

and 



Bi{x) 



3 — xi + 2x2 2j;i — X2 — 2 
2xi -X2-2 -1 + 2xi 



B2{x) 



2 + xi 1 + 3xi — X2 
1 + 3xi — X2 3 — 2xi + X2 



Then the three-dimensional spectrahedron C is represented by the LMI 

2d + 3- 



d- 

1 + 

-2 + 
-2 + 



xi + 2x2 
^ 5 + 2X2 


Sx'i — X2 
2x1 — X2 
2x1 — X2 





d + 5 + 2x2 
2d + 2 + xi 
1 + 3x1 — X2 
1 + Sx'i - X2 
-2 + 2x1 - X2 






- 2x1 
-2x1 



^2 
X2 




2d - 1 + 2x1 
d + 1 + 3x1 


1 + 3x1 — X2 





1 + 3x1 — X2 
2d + 3-xi +2x2 
d + 6 — 3xi + 3x2 



-2 + 2x1 -X2 
-2 + 2x1 --X2 

-2 + 2X1 - X2 


-; 
-; 

d + 1 + 3x1 
2d + 2 + xi 
1 + 3x1 - X2 



1 + 3x1 - X2 
1 + 3x1 — X2 
d + 6 — 3xi + 3x2 
2d + 3-2xi +X2 



-2 + 2x1 - X2 




■ 2x1 

■ 2x1 



1 + 3x1 — X2 






^2 
X2 



1 + 3x1 — X2 
2d - 1 + 2x1 
<i + 2 + X2 



-2 + 2x 



2d- 






^2 



1 + 3x1 - X2 
1 + 3x1 — X2 
d + 2 + X2 
■ 3 — 2x1 + X2 



>- 0. 



Note that the boundary of C is a surface of degree eig ht inR^ If we fix a positive 
real value for (i, then this LMI describes a two-dimensional spectrahedron. □ 
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4 Algebraic degree of the boundary 



We have seen in Theorem 13.21 that the set C is a spectrahedron. This imphes the 
property of rigid convexity, which was introduced by Helton and Vinnikov [8]. 
We briefly review this concept, starting with a discussion of real zero polynomials. 
Let f{z) be a polynomial in s variables (zi, . . . , Zg) and w a point in M*. We say 
f{z) is a real zero (RZ) polynomial with respect to w if the univariate polynomial 

9{a) ■■= f{w + a-v) 

has only real roots for any nonzero vector v £ M". Being real zero with respect 
to is a necessary condition [8] for f{z) to have a determinantal representation 

f{z) = det{Fo + ziFi + --- + ZsFs), 

where the Fi are constant symmetric matrices such that Fo+wiFi + - ■ ■ + WsFs >- 
0. Helton and Vinnikov [8] showed that the converse is true when s = 2, but a 
similar converse is not known for s > 2. For representations of convex sets as 
projections of spectrahedra in higher dimensional space we refer to [HIT]. 

A convex set F is called rigid convex if F is a connected component of the set 
€ : f{z) > 0} for some polynomial f{z) that is real zero with respect to 
some interior point of F. As an example, the unit ball {z € M" : ||^|| < 1} is a 
rigid convex set. Note that, not every convex semialgebraic set is rigid convex. 
For instance, the set {z € : zf + Z2 < 1} is not rigid convex as shown in [8]. 

The boundary of the 3-ellipse in Figured] is rigid convex of degree eight, and 
the curve (for fixed d) in Example 13.41 is also a rigid convex curve of degree eight. 
In light of Theorem 13.21 our discussion implies the following corollary. 

Corollary 4.1. The convex semialgebraic set C C R'""*"^ in is rigid convex. 

We now examine the algebraic properties of the boundary of the rigid convex 
set C, starting from the explicit LMI representation £.(x,d) >z which was con- 
structed in ()3.2p . Since the tensor product is a bilinear operation, the entries of 
the symmetric matrix C{x, d) are linear in the n + 1 variables (x, d), and they also 

(k) 

depend linearly on the entries of the constant matrices Ai and Bj . We regard 
these entries as unknown parameters, and we let K = Q(^, B) denote the field 
extension they generate over the rational numbers Q. Then K[x,d] denotes the 
ring of polynomials in the n+1 variables (x, d) with coefficients in the rational 
function field K = Q{A,B). Our object of interest is the determinant 

r{x,d) = det C{x,d). 

This is an element of K[x,d]. By a specialization of r(x, d) we mean the image of 
r{x, d) under any field homomorphism K Thus a specialization of r(x, d) is 

a polynomial in M[x, y] which arises as the determinant of the LMI representation 
()1.2p for some specific matrices Ai and B^''^ with entries in M. 

Theorem 4.2. The polynomial r{x,d) has degree NqNi ■ ■ ■ Nm, and it is irre- 
ducible as an element ofK[x,y]. Generic specializations ofr{x,d) are irreducible 
polynomials of the same degree inM[x,d]. Therefore, if the given symmetric ma- 

(k) 

trices Ai and B - are generic, then the boundary of the spectrahedron C is a 
connected component of an irreducible real hypersurface of degree NqNi ■ • • Nm- 
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Proof. As in (j3.4p we write (x), . . . , A]y^(x) for the eigenvalues of the matrix 
B]^[x). These eigenvalues are elements in the algebraic closure of the rational 
function field Q(-B, x). Here the entries of the matrix coefficient of Xi in Bk{x) are 
regarded as variables over Q. This implies that the characteristic polynomial of 
Bk{x) is irreducible over Q{B,x) and the Galois group [10] of the corresponding 

(k) (k) 

algebraic field extension Q{B,x) C Q{B,x){X\ (x), . . . , A)y^(x)) is the full sym- 
metric group S'atj, on letters. Let L be the algebraic field extension of Q(i?, x) 
generated by all eigenvalues x[''\x) for A; = 1, . . . , m and i = 1, . . . , A^'^. Since all 
scalar matrix entries are independent indeterminates over Q, we conclude that 
Q{B,x) C L is an algebraic field extension of degree Ni - ■ ■ N^, and the Galois 
group of this field extension is the product of symmetric groups S^i x • • • x -Sat^. 

The square matrix C{x, d) has NqNi ■ ■ ■ Nm rows and columns, and each entry 
is a linear polynomial in {x,d) with coefficients in Q(A, i?). The polynomial 
r(x, d) is the determinant of this matrix, so it has degree NqNi ■ ■ ■ Nm in (x, d). 
To see that it is irreducible, we argue as follows. In light of (|3.6|) . the matrices 
C{x,d) and A[L^ have the same determinant, and this is the product of the 
determinants of the NiN2 - ■ ■ Nm blocks ()3.5p of size Nq x Nq. We conclude 

rix,d) = lll[...l[det{d-Ao + xf^{x)A, + xfJix)A2 + --- + \^;2^ix)Am). 

il = lj2 = l jm = l 

The Galois group mentioned above acts transitively by permuting the factors 
in this product. No subset of the factors is left invariant under this permuta- 
tion action. This proves that r{x,d) is irreducible as univariate polynomial in 
Q{A, B, x)[d], and hence also as an (m -|- l)-variate polynomial in Q(A, B)[x, d]. 

The assertion in the third sentence follows because the property of a polyno- 
mial with parametric coefficients to be irreducible is Zariski open in the parame- 

(k) 

ters. If we specialize the entries of Ai and Bj to be random real numbers, then 
the resulting specialization of r{x,d) is an irreducible polynomial in ]R[x,d]. □ 

(k) 

Naturally, in many applications the entries of Ai and Bj will not be generic 
but they have a special structure. In those cases, the characteristic polynomial 
of Bk{x) may not be irreducible, and a more careful analysis is required in deter- 
mining the degree of the hypersurface bounding the spectrahedron C. Let hk{a) 
be the univariate polynomial of minimal degree in Q(i?,x)[Q] such that 

hk{Xram{Bk{x))) = /lfc(Amax(5fe(x))) = 0. 

holds for all points x in some open subset of M". Thus, hk{ct) is the factor of the 
characteristic polynomial det(i?fc(x) — al^i^) which is relevant for our problem. 
By modifying the argument in the previous proof, we obtain the following result. 

Theorem 4.3. The Zariski closure of the boundary of the spectrahedron C is a 
possibly reducible real hypersurface in M""*"^ . If the matrix Aq is invertible then 
the degree of the polynomial in (x, d) which defines this hypersurface equals 

m 

NQ-Wdegxee{hk) (4.1) 

k=l 
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It is important to note that, for fixed d, the degree of the resulting polynomial 
in X can be smaller than (|4.ip . Thus this number is only an upper bound for the 
algebraic degree of the boundary of the spectrahedron {x € M" : C{x,d) > 0}. 

Example 4.4. The classical ellipse consists of all points (xi, X2) € whose sum 
of distances to two given points (uj, Vi) is d. In the formulation sketched in Section 
2 and worked out in detail in [12], we obtain A'o = 1 and Ni = degree(/ifc) = 2 
for A; = 1, 2. Hence the number (|4.1|) is four, and this is indeed the degree of the 
surface c)C C M'^. Yet, the boundary of the ellipse is a curve of degree two. □ 



5 Applications and questions 

The matrix cube defined in (jl.ip has important applications in robust control [3l 
m [11] . Consider nonlinear feedback synthesis for the discrete dynamical system 



x{k + l) = F{u{k)) ■ x{k), x{k) e 
uik + l) = f{x{k))), u{k) G 



(5.1) 



where x{k) is the state and u[k) is the feedback. We assume that the matrix 
F{u) = Fq+uiFi-\-- ■ ■+UmFm G \s linear in u and that /(x) is an algebraic 

function of x = (xi, . . . ,Xm)- Such a function /(x) can be implicitly defined by 
eigenvalues of linear matrices Bk{x), that is, the coordinates of n = /(x) satisfy 

det {uk ■ - Bk{x)) =0, k = l,...,m. 

In stability analysis one wish to know what states make the system (|5.1|) stable. 
The task is to identify states x such that the feedback u satisfies ||-F(ti)||2 < 1- 
Equivalently, we need to find x such that the eigenvalues of all i?fc(x) satisfy 



Ir F{u) 



We conclude that such states x are the elements of a convex set of the form (jl.ip 
when every i?fc(x) can be chosen to be symmetric. Thus Theorem 13.21 furnishes 
an LMI representation for the set of stable states. 

In Section 2 we have seen that the LMI representation of multifocal ellipses in 
[12| arise as a special case of our construction. We next propose a natural matrix- 
theoretic generalization of this, namely, we shall define the matrix ellipsoid and 
the matrix m-ellipsoid. Recall that an m-ellipsoid is the subset of defined as 

|x G M" : ai||x - ui||2 H + am\\x — Umh < d| 

for some constants ai, 02, . . . , Om, d > and fixed foci ui, U2, • • • , Um G K""- We 
define the matrix m-ellipsoid to be the convex set 

|x G M" : Ai||x-ni||2H \- Am\\x - Unih ^ d ■ In^ (5.2) 

for some d > 0, fixed foci ui,U2, • • • ,Um G IR", and positive definite symmetric 
N-hy-N matrices Ai, . . . , Am >- 0. To express the matrix m-ellipsoid (|5.2p as an 
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instance of the spectrahedron C in (jl.ip . we fix = Ij\f and the linear matrices 



Bkix) 





Xl - UkA 











-Uk' 







Theorem 13.21 provides an LMI representation of the matrix m-elhpsoid (j5.2p . 
This generahzes the LMI given in [12^ §4.2] for the case when all Ai are scalars. 
By Theorem 14.31 the boundary of (j5.2p is a hypersurface of degree at most N2^. 

One obvious generalization of our parametric matrix cube problem is the set 



C :-- 



{x,d) 



d-Ao+E tkAk >i 0, V(ti,- 

k=l 

Amin(-Bfc(x)) <tk < Xma.AEkix)) 

for A; = 1, • • • ,771 



where Ek{x) are linear symmetric matrices different from Bk{x) for k = 1, . . . , m. 
Since the minimum eigenvalue function Amin( ■ ) is concave and the maximum 
eigenvalue function Amax ( ■ ) is convex, we see that the set C defined as above is 
also a convex set. Assuming the extra hypotheses Xmax{Ek{x)) > Amax(-Bfc(x)) 
and Amin(-Efc(x)) > Xmm{Bk{x)), the convex set C can be equivalently defined as 



{x,d) 



d-Ao+Yl tkAk h 0, V(ti,--- ,tm) : 

k=l 

XminiDkix)) <tk< Xma.ADk{x)) 

for k = 1, - ■ ■ ,m 



where Dk{x) = diag[Bk{x) , Ek{x)) is a block diagonal matrix. Therefore, (7 is a 
special case of (II. ip and Theorem 13 . 21 furnishes an LMI representation. We do not 
know whether the extra hypotheses are necessary for C to be a spectrahedron. 

An interesting research problem concerning the matrix cube (jl.ip is to find 
the smallest size of an LMI representation. This question was also raised in 
[12\ §5] for the case of m-ellipses. Theorem 13.21 gives an LMI representation 
of size NqNi ■ ■ ■ Nm- The degree of the boundary dC is given by the formula 
A^'oHfoLi degree(/ifc), by Theorem 14.31 If degree(/ifc) is smaller than A^^, then the 
size of the LMI in Theorem 13.21 exceeds the degree of dC. In this situation, is it 
possible to find an LMI for C with size smaller than NqNi ■ ■ ■ A^m? When n = 2 
and d is fixed, the projection of C into x-space is a two dimensional convex set 
Cx described by one LMI of size NqNi ■ ■ ■ Nm- The work of Helton and Vinnikov 
[8] shows that there exists an LMI for Cx having size equal to the degree of Cx- 
How can the LMI of Theorem 13.21 be transformed into such a minimal LMI? 



References 

[1] S. Basu, R. Pollack and M.-F. Roy. Algorithms in Real Algebraic Geometry. 
Springer- Verlag, Berlin, 2006. 



11 



[2] S. Boyd, L. El Ghaoui, E. Feron and V. Balakrishnan. Linear Matrix Inequal- 
ities in System and Control Theory. SIAM Studies in Applied Mathematics, 
Vol. 15, Society for Industrial and Applied Math., Philadelphia, 1994 

[3] A. Ben-Tal and A. Nemirovski. Robust convex optimization. Mathematics 
of Operations Research 23 (1998), No.4, 769-805. 

[4] A. Ben-Tal, A. Nemirovski and C. Roos. Extended matrix cube theorems 
with applications to ^-theory in control. Mathematics of Operations Re- 
search 28 (2003), 497-523. 

[5] B.-D. Chen and S. Lall. Degree bounds for polynomial verification of the 
matrix cube problem. larXiv : math . DC/0604573> 

[6] J.W. Helton and J. Nie. Semidefinite representation of convex sets. Preprint, 
2007. arXiv: 0705 .406 8, 

[7] J.W. Helton and J. Nie. Sufficient and necessary conditions for semidefinite 
representability of convex hulls and sets. Preprint, 2007. l arXiv : 0709 . 40171 

[8] J.W. Helton and V. Vinnikov. Linear matrix inequality representation of 
sets. Comm. Pure Appl. Math. 60 (2007), no. 5, 654-674. 

[9] R.A. Horn and C.R. Johnson. Topics in Matrix Analysis. Cambridge Uni- 
versity Press, 1994. 

[10] S. Lang. Algebra. Third edition. Springer- Verlag, New York, 2000. 

[11] A. Nemirovski. The Matrix Cube Problem: Approximations and Applica- 
tions. Lecture Notes, 2001, Bittp : //www2 . isye .gatech. edu/''^einirovs71 

[12] J. Nie, P. Parrilo and B. Sturmfels. Semidefinite Representation of the k- 
Ellipse. IMA Volume I46: Algorithms in Algebraic Geometry (Eds. A. Dick- 
enstein, F.-O. Schreyer, and A. Sommese), Springer, New York, 2008. 

[13] J. Nie, K. Ranestad and B. Sturmfels, The algebraic degree of semidefinite 
programming. |arXiv : math/06 1 1562i to appear in Math. Programming. 

[14] M. Ramana and A.J. Goldman. Some geometric results in semidefinite pro- 
gramming. Journal of Clobal Optimization 7 (1995) 33-50. 

[15] J.F. Sturm. SeDuMi 1.02, a MATLAB toolbox for optimization over sym- 
metric cones. Optimization Methods and Software, 11/12 (1999) 625-653. 

[16] H. Wolkowicz, R. Saigal, and L. Vandenberghe (Eds.). Handbook of Semidefi- 
nite Programming: Theory, Algorithms, and Applications, Kluwer Academic 
Publishers, 2000. 



Authors' addresses: 

Jiawang Nie, Department of Mathematics, University of California at San Diego, 
La Jolla, CA 92093, njw@math.ucsd.edu 

Bernd Sturmfels, Department of Mathematics, University of California, Berkeley, 
CA 94720,USA, bernd@math.berkeley.edu 



12 



